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Abstract 

In a hypersonic boundary layer over a wall of variable curvature, the region most sus- 
ceptible to Gortler vortices is the temperature adjustment layer over which the basic state 
temperature decreases monotonically to its free stream value (Hall & Fu (1989), Fu, Hall & 
Blackaby (1990)). Except for a special wall curvature distribution, the evolution of Gortler 
vortices trapped in the temperature adjustment layer will in general be strongly affected by 
boundary layer growth through the 0(M 3 / 2 ) curvature of the basic state, where M is the 
free stream Mach number. Only when the local wavenumber becomes as large as of order 
M 3 / 8 , do nonparallel effects become negligible in the determination of stability properties. 
In the latter case, Gortler vortices will be trapped in a thin layer of 0(e 1 ^ 2 ) thickness which 
is embedded in the temperature adjustment layer; here e is the inverse of the local wavenum- 
ber. In this paper, we first present a weakly nonlinear theory in which the initial nonlinear 
development of Gortler vortices in the neighbourhood of the neutral position is investigated 
and two coupled evolution equations are derived. From these two evolution equations we can 
determine whether the vortices are decaying or growing depending on the sign of a constant 
which is related to the wall curvature and the basic state temperature. In the latter case, it 
is found that the mean flow correction becomes as large as the basic state at distances 0(1) 
downstream of the neutral position. Next, we present a fully nonlinear theory concerning 
the further downstream development of these large-amplitude Gortler vortices. It is shown 
that the vortices spread out across the boundary layer. The upper and lower boundaries 
of the region of vortex activity are determined by a free-boundary problem involving the 
boundary layer equations. Finally, the secondary instability of the flow in the transition 
layers located at the upper and lower edges of the the region of vortex activity is considered. 
The superimposed wavy vortex perturbations are spanwise periodic travelling waves which 
are 7 t /2 radians out of phase with the fundamental. The dispersion relation is found to be 
determined by solving two coupled differential equations and it is shown that an infinite 
number of neutrally stable modes may exist. 
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Contract No. NAS1-18605 while the second author was in residence at the Institute for Computer Applica- 
tions in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. Additional 
support was provided by USAF under Grant AFOSR89-0042 and SERC. 
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1 Introduction 


This paper is the third of a series of papers reporting on our studies of the stability properties of 
hypersonic boundary layers with respect to the Gortler instability mechanism. The previous 
two papers have been devoted to the linear development of Gortler vortices in hypersonic 
boundary layers. The first, Hall and Fu (1989), is concerned with fluids which have their 
viscosity modelled by Chapman’s law; whilst the second, Fu, Hall and Blackaby (1990), studies 
Sutherland’s law fluids. In the second paper, the roles played by gas dissociation and wall 
cooling in the determination of stability properties are also clarified. In the present paper, we 
first study the nonlinear development of Gortler vortices in the neighbourhood of the neutral 
position and show how a large amplitude vortex structure can be developed under the combined 
effects of viscosity and nonlinearity. Then we consider one of the several possible types of 
secondary instabilities which the latter vortex structure may suffer, namely, the wavy type. 
For a review of the general literature on Gortler instability in hypersonic flows we refer the 
reader to our previous papers. A detailed review for the related incompressible flow problems 
can be found in Hall (1988). Here we only mention the papers which are most relevant to our 
present studies. 

The type of nonlinear theory which we use here is different from the classical weakly 
nonlinear theory and was first established by Hall (1982) and Hall and Lakin (1988) in the 
context of Gortler instability in incompressible flows. In the classical weakly nonlinear theory, 
see, for examle, Stuart (1965), the size of the disturbance is chosen such that the cumulative 
effects of nonlinearity are brought into the evolution equation as a solvability condition at the 
third order of a successive approximation procedure using a multiple scales approach. This is 
possible because the growth rate in the neighbourhood of the neutral position is small. Hall’s 
(1982b) nonlinear theory is suitable for problems which have large growth rate (correponding 
to large wavenumbers); the size of the disturbance is chosen to be so large that nonlinear effects 
come into play at the order at which the vertical structure is determined. Such large amplitude 
Gortler vortices can grow downstream of the neutral position and become so large as to produce 
a mean flow correction as large as the basic state, as is shown by Hall and Lakin (1985). When 
this happens, the vortices are confined to a core region bounded by two transition layers. In 
the core region of Gortler vortex activity, the boundary layer is forced by the vortex which 
itself is driven by the boundary layer. In the two transition layers, the vortex is reduced to zero 
exponentially. Hall and Seddougui (1989) and Seddougui and Bassom (1990) studied the wavy 
type of secondary instability which might occur in the two transition layers. Their studies were 
motivated by the experimental results of Bippes (1978) and Aihara and Koyama (1981), who 
observed that the three-dimensional breakdown of steady spanwise periodic Gortler vortices 
led to a time periodic flow with wavy vortex boundaries similar to those which occur in the 
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Taylor problem. Our present investigations are mainly aimed at finding out how the results for 
incompressible flows found by Hall (1982), Hall and Lakin (1985), Hall and Seddougui (1989) 
and Seddougui and Bassom (1990) should be modified in order to describe hypersonic flows. 

The most important property of a two dimensional hypersonic boundary layer is probably 
that it can be divided in the large Mach number limit into a wall layer and a temperature 
adjustment layer sitting at the edge of the boundary layer, and that it is the latter layer that is 
most susceptible to Gortler vortices. Since for a hypersonic boundary layer the basic state tem- 
perature varies significantly, Chapman's viscosity law is a poor approximation to the viscosity 
of the fluid and Sutherland’s law is more realistic, throughout the present paper Sutherland’s 
law will be adopted and our attention will be focussed on the temperature adjustment layer. 
It has been shown by Fu, Hall and Blackaby (1990) that when Sutherland’s law is used, the 
wall layer and the temperature adjustment layer are respectively of thickness of order M 3 / 2 
and order unity in terms of the physical variable. For Gortler vortices which have wavelength 
comparable with the boundary layer thickness (defined as the wall layer mode), the neutral 
Gortler number is a decreasing function of the local wavenumber. As the local wavenumber 
increases (physically, this may correspond to when we follow the downstream evolution of the 
vortices), the centre of vortex activity moves towards the temperature adjustment layer and 
the neutral Gortler number tends to be independent of the global wavenumber. For Gortler 
vortices trapped in the temperature adjustment layer, the neutral Gortler number is found to 
have its first term independent of the global wavenumber. This term is due to the curvature 
of the basic state; other higher order correction terms are related to viscous effects and are in 
general affected by boundary layer growth. How important the nonparallel effects are depends 
upon the wall curvature and on the size of the wavenumber. Before we review Fu, Hall and 
Blackaby ’s (1990) main findings, let us first note that for a hypersonic boundary layer de- 
scribed using a similarity variable, boundary layer growth has two lengthscales. The first scale 
is related to the similarity variable 77 (which is defined as the ratio of the Howarth-Dorodnitsyn 
variable over \plx where x is the streamwise variable). This scale, which we shall refer as the 
short scale, is not present in incompressible boundary layers and it arises because of the fact 
that r) x = 0(M 3 / 2 ) 1. The second scale is the usual one related to the variable x. It 

was shown in Fu, Hall and Blackaby (1990) that the effect of boundary growth over the short 
scale is felt mainly through an 0(M 3 / 2 ) curvature term in the y-momentum equation and that 
when the wall curvature distribution is proportional to ( 2x)~ 3/ 2 , the wall curvature exactly 
counterbalances the curvature of the basic state if the Gortler number is chosen appropriately, 
so that for this special curvature case nonparallel effects affect the stability properties in a 
similar fashion to that for incompressible boundary layers. But since the special curvature 
case is possibly of little physical relevance we shall only consider the more general curvature 
case in this paper. 
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Intuitively, the effects of boundary layer growth decrease as the wavelength becomes increas- 
ingly small, since to a very small wavelength vortex, the boundary layer streamlines would be 
almost straight lines. Therefore the effects of boundary growth are usually described in terms 
of the relative order of the vortex wavelength to the lengthscale over which the boundary layer 
growth is significant (i.e. is an 0(1) effect). Since for a hypersonic boundary layer the most 
important scale of boundary layer growth is the short scale defined above and it is related to 
the free stream Mach number, we describe nonparallel effects in terms of the relative orders of 
the wavenumber and Mach number. Fu, Hall and Blackaby (1990) showed that when the local 
wavenumber is of order unity, the downstream development of Gortler vortices is governed 
by inviscid equations and thus their spatial development has an oscillatory nature. Viscous 
effects are small but are cumulative so that they become important further downstream where 
the local wavenumber has become large. Nonparallel effects are dominant for this range of 
wavenumbers. When the wavenumber reaches 0(M x / 4 ), viscous effects become of leading or- 
der effects but nonparallel effects are still dominant. Only when the wavenumber becomes as 
large as of order M 3 / 8 , do nonparallel effects become negligible and viscous effects then dictate 
the downstream evolution properties of Gortler vortices. 

In the present investigation, our first aim is to find out how nonlinear effects compete 
with viscous effects in the evolution of Gortler vortices and to show how a large amplitude 
vortex structure can be established. We exclude nonparallel effects by assuming that the 
local wavenumber is of order M 3 / 8 . This assumption is physically relevant because a large 
local wavenumber can not only be achieved by a large global wavenumber but it can also be 
achieved by moving sufficiently far downstream. Thus a vortex of any wavenumber would 
eventually evolve into the wavenumber regime covered by the present theory as they propagate 
downstream. Our second aim is to provide an asymptotic description for the wavy type of 
secondary instability which the large amplitude Gortler vortex structure may suffer. 

The rest of this paper is divided into five sections as follows. After discussing the basic 
state and giving a brief review of the linear theory in section 2, we consider in section 3 the 
weakly nonlinear development of large amplitude Gortler vortices in a small neighbourhood of 
the neutral stability position given by the linear theory. We show that if the basic state and the 
wall curvature satisfy a certain condition, these Gortler vortices will grow until they become 
so large as to drive a mean flow correction as large as the basic state. The further downstream 
evolution of these large amplitude vortices is then studied in section 4 where we show that the 
vortices have a triple layer structure which consists of a region of vortex activity bounded by 
two transition layers over which the amplitude of the harmonic part of the vortex decays to 
zero exponentially. The position of the two transition layers are found to be governed by a free 
boundary problem which is solved in section 5 for a number of curvature distribution cases. In 
section 6 we investigate the secondary instability of the two transition layers with respect to 
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travelling waves which are tt/ 2 out of phase in the spanwise direction with the steady Gortler 
vortices. Finally in section 6 we discuss our results and draw some conclusions. 


2 Basic state and a review of the linear theory 


Consider a hypersonic boundary layer over a rigid wall of variable curvature (1/ A)k{x* / L), 
where L is a typical streamwise length scale and A is a lengthscale characterizing the radius of 
curvature of the wall. We choose a curvilinear coordinate system (x*, y*, z *) with x* measuring 
distance along the wall, y* perpendicular to the wall and z* in the spanwise direction. The 
corresponding velocity components are denoted by and density, temperature and 

viscosity by p*,T* and p* respectively. The free stream values of these quantities will be 
signified by a subscript oo. We define a curvature parameter 6 by 


6 = 


L 

A’ 


( 1 ) 


and consider the limit 6 — * 0 with the Reynolds number R defined by 

<oLP*oo 


R = 


taken to be so large that the Gortler number 




( 2 ) 


G = 2 R}! 2 6 


(3) 


is 0(1). In the following analysis, coordinates (x*,y*,z*) are scaled on (I, R~ lf3 L), 

the velocity (u*, v*, w*) is scaled on (u^, R _1/2 u^, -R _1/2 <,) and other quantities such as 
p*,T* y and p* are scaled on their free stream values with the only exception that the pressure 
p* is scaled on and the bulk viscosity X* is scaled on /£,. All dimensionless quantities will 

be denoted by the same letters without a superscript *. For an ideal gas without dissociation 
the Navier-Stokes equations are given by 





(4) 


Du dp d . du d . du 
P-RI = "ITT + 


Dt 


dx T dy^dy J ' dz^dz' 


,Dv 1_ 2 . D dp d j / \ 2 .dvp 

p( 1 — Gku 2 ) = — R— — h -7T- s (X — — p) — — 

^ Dt 2 ' dy dyV 3 dxp 


/ 2 d . dvp. 

{ 3^ dxp) + dxp ^ dy 


d . dv d dv 
+ dy^dy^ + dz^dz^ 


Dw n dp d 2 .dv 0 \ d dv 0 d dw d dw 

’■ m = - r £ + al | (A - + + P**' ’ 


(5) 


( 6 ) 

(7) 
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1 d , dT. 1 d , dT. 

?a? ( "a 7 ) + 7a7 ( '‘a7 ) ' 

(8) 

7 M 2 p = pT. 

(9) 


Here we have used a mixed notation in which (vi,U 2,V3) * s identified with (u, v,u;) and 
(xi,X2,®3) with (x,y,z). Repeated suffices (3 signify summation from 1 to 3. The constants 
7 , M and a are in turn the ratio of specific heats, the free stream Mach number and the Prandtl 


number defined by 


7 = 


M = 


ui 


U t 


00 




00 
2 1 


a == 


Moo c poo 
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where 5R is the gas constant, k is the coefficient of heat conduction, and &00 — \pr%T< 
the sound speed in the free stream. In equations (5) — (8), the operator DjDt is the material 
derivative and it has the usual expression appropriate to a rectangular coordinate system. 
The basic state is given by 


(tt, v, w ) = (u(x,y), v(x, y), 0), T = T(x, y), 


p=p(x,y), fi = fi(x,y). 


( 10 ) 


By substituting (10) into the governing equations (4)— (9) it is straightforward to obtain the 
reduced equations satisfied by the basic state. The reader is referred to the book by Stewartson 
(1964) for a detailed discussion of these basic state equations. If we define the Howarth- 
Dorodnitsyn variable y and a similarity variable 77 by 


rv y 

y= Jo pdy and 77 = Wx 

then the continuity equation is satisfied if u and v are written as 

5 " 8 “ + 

Here the functions f(r)) and T(tj) must satisfy 


( 11 ) 


( 12 ) 


//" + («*/")' = 0. (1 3 ) 

-{?$?')' + ff + /2(7 - 1 )M 2 p{f") 2 = 0, (14) 

u 

if the i-momentum and energy equations are to be satisfied. The y-momentum equation gives 


dp 

dy 


= 0 


5 


to leading order so that p = p(x). In our following analysis, we assume that there is no pressure 
gradient along the streamwise direction and therefore we can take p = constant. Equation (9) 
then gives 

pT= 1. (15) 


Once the viscosity p. is specified as a function of the temperature, equations (13) and (14) can 
then be integrated to determine the basic state. Such solutions have been given by Hall and Fu 
(1989) for Chapman’s viscosity law and by Fu, Hall and Blackaby (1990) for Sutherland’s law. 
In both cases, the boundary layer divides into two sublayers: a wall layer in which T = 0(M 2 ) 
and a temperature adjustment layer over which the temperature decreases monotonically to 
its free-stream value. Certainly, for a hypersonic boundary layer across which the temperature 
varies significantly, it is more appropriate to use Sutherland’s viscosity law 


fj.= (1 + m) 


f3/2 

(f + m) 


(16) 


where m is a constant. In this case, the thicknesses of the wall layer and the adjustment layer 
are and 0(1), respectively (in terms of the similarity variable r}). In the adjustment 

layer, the functions / and T in (12) expand as 


/ = v- 


0 , 


f(v) 

M& + \ 


T = T + 


(17) 


with / and T satisfying 


(1 + m) 


f + m 


/") +*?/" = 0, 


(1 + m) ( s/T 


T + m 


f' +7/f' = 0. 


(18) 


In the linear stability analysis, we superimpose a steady periodic stationary vortex structure 
with wavenumber a on the basic state (10) and the perturbation equations are found by 
linearizing the Navier-Stokes equations about the basic state. These linear equations have 
been fully discussed in our previous paper Fu, Hall and Blackaby (1990). It was shown there 
that for the wall mode which has wavelength comparable with the boundary layer thickness, 
the neutral Gortler number is a decreasing function of the local wavenumber. As the latter 
increases, the centre of vortex activity moves towards the temperature adjustment layer and 
the Gortler number tends to a constant which is the leading order term of the Gortler number 
expansion for the mode trapped in the temperature adjustment layer. It is this mode that 
is most susceptible to Gortler vortices since it has a smaller Gortler number than any other 
mode. 

As is typical of Gortler vortices in growing boundary layers, the evolution of Gortler vortices 
in the temperature adjustment layer is dominated by nonparallel effects. It was shown in Fu, 
Hall and Blackaby (1990) that in the hypersonic limit such nonparallel effects operate mainly 


6 


through the 0(M 3 ^ 2 ) curvature of the basic state. Thus only when the wavenumber a is as 
large as of order M 3 / 8 do nonparallel effects become negligible and the following asymptotic 
expression for the neutral Gortler number can be obtained: 


G = 


IBM 3 ! 2 


+ g 0 a 4 + a 3 


3ffo d 2 go 

. J + “ ' 1 


where 


«(*n)( 2z n ) 3/2 \/5^V 2Tq drf 

2y/2x n p^TQ 


B d =' lim M~ 3 ! 2 f°° f(r))dr), g 0 = 

M-+00 Jo 


(19) 


( 20 ) 


ctk(i„)Ti 

and To = f(r i /*), Ti = T'(^), jx 0 = £(T 0 ). The constant 77 * denotes the centre of vortex 
activity and has the numerical value of 3.001 when a = 0.72, m = 0.509. In (19) the first 
term is due to the curvature of the basic state and other terms are due to viscous effects. 
It is clear that a = 0(M 3 / 8 ) is the order at which viscous effects become comparable with 
the effects of centrifugal acceleration due to the curvature of the basic state. In the following 
analysis, our nonlinear theories will be concerned with wavenumbers in this regime. Although 
this assumption about the order of the wavenumber is necessary to obtain a relatively simple 
asymptotic analysis, it is also relevant to physically situations since our analysis is actually 
based on the assumption \/2 ~xa = 0(A/ 3 / 8 ) and the latter condition is always satisfied by 
vortices far downstream of the leading edge. 


3 Weakly nonlinear theory 

In this section, we consider the initial nonlinear evolution of large wavenumber Gortler vortices 
in the neighbourhood of the neutral position x = x n . We shall fix the Gortler number as given 
by 

^ f 2 BN \ 1 

G ~ L(* n )( 2x n )*/ 2 + g °) e*’ ( } 

where for the convenience of asymptotic analysis, we have defined an 0(1) constant N and a 
small parameter e by 

A r = M?-a“ 4 , e=l/a. (2) 

It can easily be seen that (1) is just (19) with the 0(a 3 ) term on the right hand side neglected. 
This implies that the linear neutral position corresponding to (1) would be x n + O(e), In the 
present context, it does not matter whether we start with (19) or with (1). 

We now turn to the derivation of the nonlinear perturbation equations. The total flow is 
written as 

u~u-\- — U } v = v + F, w = w, 

Mi 

P = P+j(p+P)> T = T + T, (3) 
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where Mi defined by 
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Mi = — rrr 

M5? + i - - 

is used to scale the streamwise perturbation velocity so that M will not appear in the following 
nonlinear perturbation equations. In (3d), p/ R is the second order correction to the basic state 
pressure (the leading order term p is a constant). On substituting (3) into the Navier-Stokes 
equations (4)-(9), making use of (12) and (17) and neglecting cubic and higher order terms 
of the perturbation quantities U,V,W,P and T, we obtain the following set of perturbation 
equations: 


l_.dU_ _ V_dU_. 
r dx 2x dr) 


pU zz - 


9 -d B -^) + 


f" 


2xfdrj K fdr) y/^xT 2 
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2pT' ... U d 
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1 dV 
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T dV rj_dV 
f 2 dx 2x dr) 


4p dV 




dT 


-I 
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4A 
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y/2xT dr) 
p dW z 


3 • (2x)T 2 dr) dr) 3 ■ (2x)T drj K f drj J 3 y/2xfdr) 
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1 dW 
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Here p = p(T )//. 2 = dpjdT,p = d 2 p/df 2 and the bulk viscosity has been taken to be zero. 
The linear perturbation equations discussed in Fu, Hall and Blackaby (1990) can be obtained 
from the above equations by neglecting nonlinear terms. 

To solve these nonlinear perturbation equations, we first note that in the large wavenumber 
limit, Gortler vortices are trapped in an internal viscous layer of 0(e 1 ^ 2 ) thickness centred at 
the most unstable position 77 = 77 *. We therefore define a variable <f> by 


4> = 


y-v 

C l/2 


(9) 


In the neighbourhood of the neutral position x = x n> the growth rate in the x-direction can 
be shown to be of order 1 /e and it is appropriate to describe such rapid growth by defining 
another variable x by 

x = ( 10 ) 


An order of magnitude analysis of the perturbation equations (4) — ( 8 ) then shows that 

U = 0(e 2 V), W = 0(e 1/2 V), T = 0(e 2 V), P = 0(6~^ 2 V). ( 11 ) 
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We assume that nonlinear effects reinforce the fundamental at the order in e at which the 
vertical structure is determined. It can then be deduced that the appropriate sizes of the 
perturbation quantities must be 

U — 0(c 3 / 2 ), y = 0(e" 1 / 2 ), W = 0( 1) 

T = (e 3/2 ), P = 0(c' 1 ), (12) 

and that the corresponding mean flow corrections, signified by a subscript m , must be of the 
orders 

U m = 0(e 3 / 2 ), V m = 0{e), P m = 0(e- 2 ), T m = 0{e 3 ' 2 ). (13) 

As is well known, the nonlinear interactions that occur in the Taylor-Gortler problem do not 
generate a mean flow in the spanwise direction. We therefore look for asymptotic solutions of 

the form 

U = e 3/2 |u to0 + c 1/2 Um i + H [(Uo + e 1/2 Ui + • • -)E -\ h C.C]} , 

V = e- 1 ' 2 {(Vb + € 1/2 Vx + • • -)E + ■ • • + C.C .} , 

W = (Wo + e 1/2 Wi + ...)£ + ••■ + C.C., (14) 

T - e 3/2 {0 mO + e 1/2 0 m i + • • • + [(*o + c 1/2 ^i + • • -)E + • • • + C.C.]} , 

P = e- 2 (P m0 + e 1/2 Pm i + •••) + * _1 {(^o + c 1/2 Pi + •••)£+••• + C.C.} , 


where 


t lz \ 

E = exp(— ), 


(15) 


and C.C. denotes the conjugate. On substituting these expansions into the perturbation equa- 
tions (4)-(8), expanding all coefficients there about x — x n and T) = Tj*, and then equating the 
coefficients of like powers such as f, Ee and etc., we obtain an infinite hierarchy of equations. 
To leading order, the Gortler number g Q is determined from a solvability condition for (Vo, Po) 

and is given by 

2 v /2i^Ao 7 o 


go = -■ 


CTK 0 i i 


(16) 


whilst Uo,Wo,6 0 and Po are related to Vo by 

fs 


U 0 = - 


iW o = - 


^/2x n fioT 2 
1 3 Vo 


zrVb, 6 0 = - 


l 




V2^T 0 d* ’ 


Po- - 


V^^T 2 

po 3 Vq 


v^^o d4> ' 


(17) 
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where fg = T 0 = T(rf), T x = f'(r, f), Ao = A(^o) «o = *(*„)• Note that (16) is 

consistent with (20), as we would expect. The mean pressure, the mean stream wise velocity 
and the mean temperature are found to satisfy 


9P mO / W2l n g 0 

~2T 0 


&m0 = 


(18) 
(19) 

PZTl H ' (20) 

To next order, we obtain four expressions similar to (17) for Z7i , 0i, W\ and Pi in terms of 
V\ and Vo and the condition that 

( 21 ) 


and 



dcf> 


2z n <7T 0 

d 

d 2 

k Ao 

dx 

d<t> 2 

/ 2x n To 

d 

d 2 ’ 

V AO 

dx 

d<j> 2 


^mO — 


^mO) 

2<j 2 Ti a|vbi 2 

rtf 2 d<t> ’ 

2 fg d\v 0 \ 2 


d 9o, - n 

■^1,=,- - 0, 


which implies that rf is where g 0 attains its minimum (and hence that rf is the most unstable 


position). 

If we carry out the expansion to one order higher, we find from a solvability condition for 
(V 2 ,02) that V 0 must satisfy the evolution equation 


d 2 V 0 

d<t> 2 


2(1 -f < 7 )fpx n dVp 
3Ao dx 


— a<j) 2 Vo + bxVo — 


2x n T 2 Tr dQmo 
37i ° d<f> ’ 


( 22 ) 


where 


a 


Tp z n d 2 g 0 
Zgo dr] 2 


> 0 , 


- 2z n Tp f 2 BN_ . K\ 3 . «i 1 \ /23) 

3 1 5 o«o( 2 z n ) 3 / 2 ' Ko 2z„‘' Ko 2x n J ’ 

and where K\ = Thus we have two coupled evolution equations (19) and (22) together 

with an uncoupled equation (20) which govern the downstream evolution of the fundamentals 
and the mean flow corrections. It is easily seen that the present weakly nonlinear theory is 
different from the classical weakly nonlinear theory (see, for example, Stuart (1965)). In the 
latter theory the amplitudes of the mean flow corrections are one order smaller than those 
of the fundamentals and thus the resulting evolution equation for the fundamentals is an 
ordinary differential equation and the former can be determined independently of the mean 


flow corrections. 

Before we present the solutions of these nonlinear evolution equations, let us first note that 
equation (22) reduces to the linear evolution equation 


d 2 v 0 

d<t> 2 


2(1 + a)foXn dVo 
3Ao 9x 


— a(p 2 Vo + bxVo = 0, 


(24) 
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after the nonlinear forcing term on the right hand side has been neglected. This equation has 
been discussed in Fu, Hall and Blackaby (1990) and its general solution can be obtained by 
first looking for separable solutions and then expanding in terms of the eigenfunctions. It is 
shown there that the mth mode is neutrally stable at x n = (2m + 1 )y/d,/b. The most unstable 
mode corresponds to m = 0. Therefore Gortler vortices with Gortler number given by (1) are 
neutrally stable at x n + ex n . If we replace x n in (1) by x n — ex n and expand the resulting 
expression up to and include the 0(e~ 3 ) term, we recover (19) which is the appropriate Gortler 
number expansion for Gortler vortices neutrally stable at x = x n . Thus, as we remarked in the 
paragraph below (2), it does not matter whether we use (19) or (1) for our weakly nonlinear 
theory; such a difference in the choice of the Gortler number only results in an 0(c) shift in 
the linear neutral position. 

We now discuss the solutions of the nonlinear evolution equations (19) and (22). We shall 
not consider (20) any further since its solution for the mean streamwise velocity u m o is not 
needed in the remaining discussions of this paper. To simplify the notation, we make the 
following substitutions: 

X = iU C = (4a) 1 / V, 


V 0 = - 


# 2x 

/i 0 (4a) 1 / 4 V *~3 


2Va 


Lo = — 


2s»f 2 

6T!(4a)i/4 


mO • 


Equations (19) and (22) then become 

( d 2 2d 


d( 2 1 + a dX 


&m0 = 


5|Vb| 2 


d( 


<L-hJL-lt 2 + X 

d( 2 3 k dX 4^ ± 


where 


k = 


) Vo = 
(i + g)x n ro|fe| 


oi> d®rn.O 

2Vo ~dT' 


(25) 


(26) 


(27) 




In (27), the positive (negative) sign is to be taken if b is positive (negative). These two 
equations are of the same form as Hall’s (1982b) equations (3.15a,b) for incompressible flows 
(they are identical when a = 1). The reader is referred to that paper for a detailed discussion 
of their numerical and asymptotic solutions. According to Hall (1982), an important property 
of these two coupled evolution equations is that any initial disturbance introduced upstream 
would either decay to zero or evolve into a unique large amplitude structure at large downstream 
locations, depending on whether b < 0 or b > 0. The latter conditions are in fact the conditions 
for Gortler vortices to decay (6 < 0) or grow (6 > 0) linearly downstream of the neutral position 
+ € %n (see Fu, Hall and Blackaby (1990)). 
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The easiest way to see how a large amplitude structure is possible is by looking for such 
solutions for (26) and (27) directly. If an asymptotic state is to be achieved, it must be the 
last two linear terms on the left hand side of (27) that balance with the nonlinear term on the 
right hand side of (27). Thus after integration we have 

Marsel (28) 

where, based on the argument given in Hall (1982), we have assumed 0 m o to be an odd function 
of £ so that we can put the arbitrary integration constant to zero. On substituting (28) into 
(26) and integrating, we obtain 

{*-<£>■}• (29) 


where C is an integration constant to be determined. Solutions (28) and (29) are only the first 
order approximations. We note that the similarity variable £/ is important. Thus if we 
define 


£ 


\/~X 


and look for the following form of asymptotic solutions for (26) and (27): 


(30) 


6 m 0 = -^ 3 ^ 2 ^mOo(0 + X l t 2 QmO\{£) + ' ' ' » 


Vo = x 1/a vbo(0 + *~ 1/2 Voi(0 + • ■ • , 


(31) 


we would have 6 m oo and Voo given by 


(32) 


Here, without loss of generality, we have assumed Vo to be real. 

The above solution breaks down near £ = ±<7 since Vqq must be positive. In each of 
these regions Vbo develops a boundary layer structure (hereafter we shall call them transition 
layers) and it can be shown that each of the required layer is of thickness X -1 / 6 . In the upper 
transition layer £ = C we define a new variable 0 by 


0 = X l/6 (CX l/2 - C) = x 2 t\c - o* 


(33) 


To match with the core region solutions (31) and (32), we have to look for solutions of the 
form 


0mO = + * 5/6 *mo(V0 + + ’» 


Vo = X 1/6 Vo°(V0 + X~ 1/2 V C f(V0 + • • • • (34) 
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On substituting (34) into (26) and (27), equating the coefficients of like powers of X and then 
considering the matching conditions at various orders, we find that and are simply 
the expansions of the core region solution (32a) in the transition layer, i.e., 


and that Vq 5 satisfies 




7T = c <^ + tt^° - 2D *° ~ 2(<w - 


By matching with (32b), we require that 

(Vo °) 2 




(35) 


(36) 


(37) 


Equation (36) is a particular form of the second Painleve transcendent. The existence and 
uniqueness of its solution has been proved by Hastings and Mcleod (1980), and this solution 
is proportional to when ^ -oo. Thus % decays to zero exponentially in the two 

transition layers which bound the core region of vortex activity. 

Above the upper transition layer and below the lower transition layer, (26) reduces to 




It admits an asymptotic solution of the form 

0mO = X 3/2 0mC>(0 + 


To match with (34a) and (35), 0 mO (O must satisfy 

dOmO 


0mO 

Substituting (39) into (38) gives 


i(C-±C% 


d( 


1 

2 


91 

8 




( d 2 ak d 3cr \ ... 

[d? + T+7 ' f d? ” T+cr k ) mo(0 " 


(38) 


(39) 


(40) 


(41) 


The solution of this equation which decays to zero when ( — * oo is 

0 m o = -Aexp I - I • U 


( 4 ) 


7 

2’ 



(42) 


where 17(7/2, yjkj 2£) is a parabolic cylinder function, A is a constant and k is defined by 

2 a 


k = 


1 + cr 


-k. 


(43) 
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On using (42) in (40), we obtain 



1 - C 2 / 4 
C-C 3 / 12’ 


A = 


Cj-c^/n 

2Z7(7/2, yfk/2C) 


• exp 



( 44 ) 


The first of these relations determines C, the location of the upper transition layer. It is easy 
to see that if we multiply both sides by the denominator of the right hand side and then move 
the numerator to the left hand side, the resulting algebraic equation changes sign at places 
where the denominator and the numerator vanish; so there is a single solution which lies in 
the interval 2 < C < 2\/3 for all k. 

The solution for the lower transition layer can be obtained in a similar fashion. It can be 
shown to be given by (42) with A replaced by —A and ( replaced by Thus by (31), (32), 
(34), (37) and (42), for a given X > 1 the flow structure in the interval — oo < 4> < oo iB 
completely determined. 

Since the amplitude of Vo and 0 grows as the vortices propagate downstream, insertion of 
(31) and (32) back into (14) shows that the latter expansions become invalid when x = 0(e~ l ) 
where the mean flow corrections € 3 ^ 2 u m o and e 3 ^ 2 d m o become as large as the basic state. Since 
<f> — — 77*), at x = 0(e“ 1 ) (i.e. x — x n = 0(1)) the transition layers are at 77- rf = 0(1) 

and are of thickness of order e 1 ^ 2 ^” 1 / 6 = e 2 / 3 . In the next section, we shall consider the further 
downstream development of these large amplitude vortices beyond x — x n = 0(1). 


4 The fully nonlinear theory 

It has been shown in the previous section that at positions 0(1) downstream of the neutral 
position x n > the mean temperature correction and the mean stream wise velocity become as 
large as the basic state. When this happens, we expect that the large M structure of the 
boundary layer is still valid. The total flow is now written as 

u — u + -tt-U, v — v + V, w = W, 

Mi 

P = P+j(P+P), T = T + T. (1) 

It should be noted, however, that although (1) is of the same form as (3), (u, v, p, T) here 
is the non-harmonic part of the total flow and is different from its counterpart in (3) which 
represents the unperturbed basic state. In the temperature adjustment layer, the similarity 
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variable 77 = rj(x,y) is defined by 


1 P 

V ~ yfeJo 


dy 


( 2 ) 


T(x, r)(x,y)y 

where the function T in the integrand is understood to be the composite solution of the mean 
temperature (i.e. the wall layer temperature plus the mean temperature in the temperature 
adjustment layer). We note that because of the 6(l) correction from nonlinear interaction, the 
mean temperature is now also a function of x. In the limit rj — ► 0 or x — x n — » 0, T(x, rj) — * T(rj) 
and (2) then reduces to (11). 

We assume that the mean velocity components u and v have the following expressions: 

. = a/M 5= _i={-f(x, ,)/+?!/(?)} + »,(«,-,). (3) 

Here /( x,tj) and f(x,-q) expand as 

P 


+ %^+ 


f = f (*,!/) + 


(4) 


f( x ,v) = v- M \/2 ' Mi 

The function v s (x,r)) in (3b) is added in order to satisfy the continuity equation. As can be 
seen from (31), this added term is partly due to the dependence of T on i and partly due to 
the 0(1) me an flow correction from nonlinar terms in the continuity equation. The function 
I(T) in (3b) denotes the integration of the mean temperature from 0 to 77. It includes the 
contribution from the integration of the wall layer temperature and thus has the expression 


J(T) = m^ 2 b + £ (f(x,0 - d Z 


3(1 + m) 2 


er 2 r ) 3 


(5) 


where the constant B is defined by (20a). 

With the aid of (2)-(4), the following important relations can easily be established: 

1(f) 1 dl(f) 


Vx = ~ 


d _ d 

*te + V di 


2 xT 

f v s 


T dx 


_ d_ I _vi_ _j__ <L 

dx + \</2 xT 2x T dx j drj’ 


1 df 

U ~ 1 + Mi dr) + 


V = 


*/2x 


M 3 ' 2 + D(x, 77 ), 


where D(x,r 7) is defined by 
1 


D <*■ = vb {f { f{x '° “ + ^ " 


vf 


( 6 ) 

(7) 

( 8 ) 

(9) 


V2i\J 0 v*'”' 7 J V®' 

The operator on the left hand side of (7) will frequently appear in our following analysis. To 
simplify notation, we shall denote it by L(). Thus for any function F(x,rf) we have 

_ dF_ f vs V 1 dI(T) \d£ 


L(F) 


dx + \ y/2xT 2x T dx j dr}' 


( 10 ) 
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By substituting (1) into the Navier-Stokes equations (4)— (9), making use of the above 
relations and then neglecting all the cubic and higher order nonlinear terms of ( U^V^W^P^T ), 
we obtain the following perturbation equations: 


1 L (r\~ J l-JL(!£-)-1 L (v}-uu l - d ( i t du \\ v 

T L "> 2xTd 1 / T ' T™ ' * “ 2xTdrfT dt)' + y/iif 2 

+ srkirw - + ? r % + 

-KTu.).-T l m + -rLv ] - ^4^ 


T 3 AxTdr)'' T 

dp 


T Uri> + 2T ' G ' 1 “ 1 + V^If Stj 6if arj' T 


4 *(£) 


1 nT f 4 <9 7# 

+ y £ (V) - 2x f2 F _ " 3 (2 x)fl$Tlhp + 

f 1 , 5 ,,, 2T d T 2 llv' dT 

| -GKU - ^3 J2 M + ( ) + 37^( -f)j f2 ~ 3x f2 Q n 


1 dP 
\/2xf dr] 

2jiv’ dT 


2uT' 

+ -4 =^W Z - 


3\/2xf 3v/2 xT dr) 


d W*) + hnL*y%: + wv a ) 


T K sf2iT dr\ 


v' 


T 2 ( L (V) + 

4 dJidV 


_ - v r 2 . 4A flv. i ar 

+ [ 3 MWz 3v^r dr ) 1 ‘ V2if dr) 


-{ 


+ 


4£ 


\ m a_ / T 


5W, 


3 • (2x)T 2 dr) dr) 3 ■ (2x)T dr) K T dr) J 3^2 xfdr) 


dp 


3\f2xTdr) 
r£ 2 


w. + AVz 


! } T [ V2i 


dW 


+ + ^Gku 2 ] 


l V2Zf dr) + VM 

_LAfKr 2> i =0 

6xTdr)\2T ) 


pT 


^V z + 


%f2 xT z 3-<j2xT dr) 


dV L _ p , t _ 1 ( dW_±dW 

’ o t C\ \rrr^ z /tT\ 


3 • (2x)T 


—V 


T v dz 


4 jj z i d ,jidw s I, 

+ 3 Wzz + teT dr, [ T dr) > T [ '^f’dr, 


T_dW 

rf>2 V 


7] dw 


T 2K dx 2x dr) 


) + U 


pW zz + 


dv z + 1 


2x dr) 

+ WW Z ) 


) 


3y/2xT dr) *J2 xT dr) 


idpdw j^^idw 

' ok m 1 ) O O * m Cl \ m r\ J f ’ 


2xT 2 dr) drj 2xTdr) y T dr) 


■>} 


W 9T 

v/alr'VSr 3>j ‘ ; 3>i 


*,(■ 1 
rrr V 


4 2 

+ “ 3V2xf dr) 


9V. mm n 

)pT z = 0, 


T* L ^ + V^f 2 + V2Zfdr, { T )+ T Wz s/2it*' T ^ T 2 ) 


( 11 ) 


( 12 ) 


( 13 ) 
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( 14 ) 


_ f -I—™ +W)^~ l-L-V*- + Wri-vl 1) 

( V2Jf dr, +Wz) f* dT, + VV dz A T* ) 


T 2 v? 

+ L W + v ^ T=0 ’ 

d _ _ r _ 1 _ f £(£) , 

?7 7 i^ r ) + v / 2^r2 F+ T jL( } 1 T2 + 


2xaTdr\ \T ) 
P-T' dT 


2xoTdr) K T n 


-Ult — 1 , — (jfl— ) [ 1 ( + wr,) 

<7 “ 2xT 2 o drj 2xaTdr, K T drj’ T K s/2if dr) 

T T' 1 dii dT 

~ 2^T 2 d^ T ^i 

_ Ar J, ±/T AT _ , , HQt 2 - JL, A ( EIt 2 ) = 0. (15) 

a l SfAj l f^ J + lii,W+ T» 4 <tiT dr)\ T 1 ) K) 

Here, to simplify notation, we have used a prime to denote partial differentiation of the mean 
flow quantities with respect to rj . As x - x n -+ 0, the mean corrections produced by nonlinear 
interaction become increasingly small and the above equations then reduce to two sets of 
equations: the basic state equations (13) and (14) and the perturbation equations (4)— (8). 

In the light of the results given in the previous section, we expect that Gortler vortices 
would be trapped in an 0(1) region bounded by two transition layerB centred at tj ~ tji (x) and 
7 ft(x), each of which has thickness of order 0(e 3 / 2 ). The configuration is sketched in Figure 1 
in which the region of vortex activity is denoted by I, the upper and the lower transition layers 
by Ha and Ha respectively, whilst the region above the upper transition layer and the region 
below the lower transition layer are denoted by Ilia and Illb, respectively. The flow properties 
in these regions are now considered separately. 

We start with the core region L There the sizes of the perturbation quantities can be 
determined from the results given in the previous section. From (25) and (31) we deduce that 
at x — x n = 0(1), 


Vo = 0(e- 1 ' 2 ), ^ = 0(1). 


Relations (17) then give 


U 0 = O(e- 1 / 2 ), B 0 = 0(e- 1/2 ), W 0 = 0(1), P 0 = 0(1). (16) 

Then from the relations between (Uq, V»> Wo, Bo, Po) an< l {U> V, W, 6, P ) shown in (14), 
we deduce that 

U = 0(e), V = 0{e~ 1 ), W = (1), T = 0(e), P = O^' 1 ). (17) 

We therefore assume the following form of solutions for (11)— (15): 


j- = fo(x,rj) + ef 1 {x,r)) + ---, 
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( 18 ) 


T(x,r]) = To(x,y) + H > 

v s = v$(x,ri ) + »}(*, + 

If = c{e(^ +eU} + ••■) + iE 2 {Ul + •••) + ••• + CC.}, 

V = c _1 {p^ 1 + eVj 1 + ••■) + +•••) + "• + C '- C7 -} > 

w = {e(Wo + «w? + •••) + *e 2 (wI + •••) + ••• + CC ) ’ 

P = c - 1 {E(Po 1 + 6P 1 1 + -*0 + ^Vo+-") + --+ C ' C, }» 

T = e {-E(^o + e ^i + "*) + f-E 2 (^o + ■*’)+ ^ C.C.} , 

where £ is defined as in (15) and C.C. denotes the conjugate. We now substitute these 
expansions into the perturbation equations (11)-(15). After equating the coefficients of E e 
in (14), (11) and (15), we obtain 

dfo v° s dfo 

+ s/2xTq dr) To dr) dx yfe dr) 


2T 0 d ., rl . 2 


\n xH{x)dr] 


(AolVb 1 ! 2 ), 


(19) 


L(*L d Jl\ + JJ°- 2«^+(2x 

dr) \T 0 drj J V dr) dx \ 


dI &l -VTxv's 


dx 


7 T 


dfo 


To dr) 


+ 2xU ° ■ + ^ 


adr)\Todv) V dr) dx \ 


y/2x dfo i t7 i 

T 

a/ (^o) _ o\ 

VZl 4 J To dr) 


(20) 


dx 


= ^ V S f + 2«*J • + cc - 


( 21 ) 


where pc = fi(To)> and a bar over U and W signifies conjugation. 

On equating the coefficient, of Be' 1 in (11), Be" 1 in (12), Be" 1 in (13), Be’ m (14) and 

Ee- 1 in (15), we have 

^ + 7m^ v °' =0 ’ 

_ T/ 1 ^(*) 0 1_ O 

W)Vo m2~°0 ~~ U> 

J 0 

+ IM^o 1 ) = 0, 

72if 0 dr) 0 T 372130 dr) 3 


A(S) + vS ( iw„') = 0 , 

1 07o T _i 1 _ /)i _ a 

- 7 — -r - 5 — Vo + “M> — 0 > 

\/2iTq drj a 


( 22 ) 

(23) 

(24) 

(25) 

(26) 
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where Jiq = /z(To) and where H(x ) is defined as 

H{x) = BN { - ——p J + - K g 0 . (27) 

Equations (23) and (26) are consistent only if 

°_STo = „ (28) 

jSplo H M 

This is a first order differential equation which can be solved to give an expression for the mean 
temperature. Note that in the limit x — > x n , (28) reduces to (16) which is the condition for the 
original basic state to be neutrally stable at a single point rj = rf . Equation (28) shows that at 
0(1) distance downstream of the neutral position, the basic state is forced by the vortices to 
be such that it is neutrally stable everywhere simultaneously in the region of vortex activity. 
With the aid of the relations (22)-(26), we can express Ul,Wl,8\ and in terms of V£ 


as follows: 


el = - 


\Z2x{1oTq dr} 


dT °V l U 1 - 1 d/o T ^1 

°’ °~ V^toTSdr, °’ 


iwr » - f >■ 

F ° = + + wj 

With the use of these relations, equation (19) reduces to 


(30) 

Integrating this equation then gives 

r? = v^^|^) + f^|V 0 '| ! . (31) 

Here we have put the arbitrary integration function of i to zero. Even if we had not done so, 
this function could be shown to be zero at a later stage. 

With the aid of the relations (29)-(31), we can reduce (20) and (21) to 


d (nodf 0 \ df 0 dfp 2 d ( 1 d fo_ lv i,2 

fo, 1 % d V I + V d V dx ~ Tq dr) y JlqTq dr) 


4\/2x ^/o| T ^-i|2 ( 

7 h { x ) tk> ^ v ° 1 ’ (32) 


1 d f}iodfo\ dtp dfp 
a dr) l T 0 dr) J ^ dr) * dx 




After To(x , 77) has been determined from (28), equation (33) can be used to determine |Vq | and 
solving (32) then gives an expression for /o(x,7/). From (28) we see that the boundary layer 
flow is forced by the vortices which, from (32), are driven by the boundary layer. 
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given by Sutherland’s law (16), (28) becomes 

dfp = 
dr] 


rp7 

1 0 


Integrating (34) gives 


1 ( l m * I K* I 1 V*\- ( 1 + m ) 2 

+ 5 mTo+ 4 T o) - “ 


Tq 6 


where a(x) is a function to be determined. 
Integrating (33) with the aid of (34) yields 


IW = 


H{a 


l temperature To(x y r]). 

When /x is 

V2x 

H(xY 


(34) 

JTx ( ^ 

H(x) {T,+ 

a(x)), 

(35) 

df 0 _ df 0 \. 
■’ih, 2x -^j dri 

(36) 


2\/2 xfioTo J To 1 

We now assume that |Vq | 2 vanishes at 77 = 771 ( 1 ), 772 ( 1 ), which bound the region of vortex 
activity. As in the case of weakly nonlinear development discussed in the previous section, 
these two boundaries are also where the solution (36) breaks down and where transition layers 
exist. The thickness of each of these layers is 0(e 2/3 ) so that in the upper transition layer at 
77 = 772 ( 1 ), we define 




V~ V2 

,2/3 * 


(37) 


Near 77 = 772 , we can deduce from (36) and (29) that Vq , and Uq are all of order e 1 ^ 3 , and that 
Wq and Pq are both of order e -1 / 3 . With the use of these results, the sizes of ( U,V,W,T,P ) 
can be determined from (18), which are shown in the following asymptotic expansions: 


u -- 

= e 4 / 3 { 

E(U 0 i 

+ e 2/3 U u + • 


• + C.<7.} 

V = 

- £ -2/3 

{^(Voi 

+ e 2/3 V n + ■ 

•*) + ■■ 

■ • + c.c) 

w = 

6- 1 / 3 ] 

[E(Woi 

t + c 2/3 W u + 

•••) + 

• ■ • + c.c. 

p = 

€“ 4 / 3 • 

{e(p 01 

+ e 2 / 3 Pu + • 

••) + •• 

■ + c.c) 

T 

= £ 4 / 3 j 

[E{6 01 

+ e 2 / 3 0 n + • • 


• + c.c) 


(38) 


^ = fo(x,Q + e*'*h{x,0 + e 4/3 / 2 (*,0 + • • ■ , 

r = f 0 (*,0 + c 2 / 3 ri(x ,0 + c 4 / 3 r 2(*,0 + ---. 

Vg = vso(x,£) + e 2 / 3 i«i(x,£) + e 4/3 vs 2 (x,() H . 

Here we expect that the first two terms in (38f, g, h) are simply the expansions of the mean 
flow functions f 0 (x y 7]) } To(x y 7]) and in I near rj = 772 . Hence 

fo(x,0 = M x >V 2 ), h( x >0 = g ^°^ T? ,T?2 ~ 
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(39) 


f 0 (*,0 = fo(*,ift) T 1 (x,0 = — 

Similar expressions can be written down for £$o and fyi« On substituting (38a-g) into (11)— (15) 
and equating the coefficients of various orders of the harmonics and mean flow quantities, we 
obtain a hierarchy of equations. From (12) we have from equating the coefficients of Ee ~ s / 3 
and Ec~ 2 , 

B(x) I 


MO^Ol 


n 


-0oi = 0, 


(40) 


Mo^ii + £o£7oVoi jjr^oi) 


4/zo d 2 Vo\ 


+ 


n 

dPoi 


To 

no d(iWoi) 


= 0. 
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Here and subsequently we write Tq for dfo{x, T) 2 )/dr) to simplify the notation and Jiq = 
p.(To(x, 772)), Jio = KTo(x,Tf 2 )). Equating the coefficients of Ee~ 7 ? 3 and of Ee~*/ 3 in (14) 
gives 

T. _ j-i 1 . A 

(42) 


Mo 


z^Ji ? W- Pm + t MiWm) = 0 ' 


dv °\ + = o 


s/2xT 2 d( ' T 0 

Finally, by equating the coefficients of Ee~ 2 ^ 3 , Ee° and E°e° in (15), we obtain 


(43) 


rpf 
1 0 

V2xT 2 


Voi + —0oi = 

a 


0, 


n 


Vu + + 75 ~ ^t )Vc 

+ vW°'-^ 9 -w- = °- 


01 


V2T0 _ mo ^ 2 r 2 _ 1 / - _ Ao\-f/2 

2iT 0 2xcrT 2 d( 2 2x oT^ T 0 ’ 0 


(44) 


(45) 


+ ^^(F 01 ^ + Voi ^1) + ^(^oi0oi - i^oi0oi) = 0. (46) 

We have not written down the equations obtained from the x- momentum equation since the 
determination of Vbi , TVoi , 0oi and P 0 i does not involve the streamwise velocity component. 
Equations (40)-(46) can be solved in the following way. First, from (42), (43) and (44), we 
can express iWpi,Poi and #oi In terms of Voi and its derivatives. We note that equations (40) 
and (44) are already consistent because of (28). Next, equations (41) and (45) give a 2 x 2 
inhomogeneous matrix equation of the form A( = f for ( = (Vu 7 9u ) T . The inhomogeneous 
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term f involves dT 2 /d£ as well as Voi and its derivatives. The former can be determined from 
(46) and can be shown to be given by 


df 2 (jTqt^T'q , Ao 1 2 < 72 ^|T' r n, I 2 - S(x) 

W~ T 0 ] °^ 5()> 


(47) 


where S{ i) is a function to be determined. Since the coefficient matrix A has zero determinant 
because (40) and (44) must have non-trivial solutions, the inner product of the left eigenvector 
of A and f must vanish. Omitting all the details, we can show that the latter condition can 
be reduced to the following second order partial differential equation for Vqi : 


d 2 V c 


01 


4x 


d ( 2 


+ Si(x)£F 0 i = — • ry VoilFoil 2 + S 2 (x)V 0U 

o fl 0 


where 


Si(x) = -^f 0 (*?' + f V2T 2 ) , 

O \ flQ 1*0 / 


(2s) 

3 


*(*) = ^ ■ £ ’ Si*)- 

1 0 


(48) 


(49) 


(50) 


This equation is a particular form of the second Painleve transcendent and has been shown by 
Hastings and Mcleod (1978) to have a solution such that 

4xa 2 


3^o 


■|Vbi| 2 ~ Si(x)£ as £ — ♦ -00 


(51) 


and that |Voi| decays to zero exponentially as £ — ► oo. Thus the Gortler vortices are trapped 
below region Ilia and the condition (51) ensures that V 0 \ matches with VJ in the core region 
I. An identical analysis applied to the lower transition layer lib shows that the vortices there 
are also reduced to zero exponentially away from the core region so that they are also trapped 
above region IIIb. As a consequence, above the upper transition layer and below the lower 
transition layer, there are only mean flow fields. 

In IIIa,b the mean flow fields are still formally represented by (3) and (4), but now / and 
T expand as 

f(x,v) = fo(x,v) + 0(e), T(x y r}) = f 0 (x,T)) + 0(e), (52) 

where fo(x,ri) and To(i,t?) satisfy 


„ , / y/fo d 2 f 0 \ , J 2 fo 0 . d 2 f 0 

(1 + m) dj [fo + mdr, 2 ) +T1 dTr> 2x dxd v 


= o, 


(1 + m) d I yTo dTo 


dy\f 0 + mdrij ^ dt) * dx 
To satisfy the continuity equation, vg must be calculated from 


vs 


= \plx 


dI(T) 
dx ’ 


(53) 


(54) 
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where the function I is defined by ( 5 ). We note that the governing equation ( 53 b) for f 0 (x,r)) 
is decoupled from that for /o . We can therefore solve ( 53 b) on [0,771], [772,00] first subject to 
the following boundary and matching conditions: 

(0- 

* 9 1 

T 0 (x, 77) — » — =■ (1 + tto) 2 -7 + ■ • • , as 77 — ► 0 , 20(1,77) -» 1 , as 77 — > 00; ( 55 ) 

a 2 rj* 

(ii) . 

To(x,Vj) = To{x, rjj), df 0 (x,T}j)/dTj = dT 0 (x 1 ]j)/dr], j = 1 , 2 ; ( 56 ) 

(iii) . 77 = 771 and 772 are where Vq vanishes and from ( 36 ) they satisfy 


/■» 1 Jl d.podTo df 0 dTol 

L T~o \<Td V { T 0 d V ) + V dr ] lX dx ] 


drj = 0. 


( 57 ) 


Thus equations ( 52 b)— ( 57 ) constitute a free boundary problem which can be solved numeri- 
cally for a given curvature distribution to determine the two boundaries 771(2), 7fc(x) and the 
unknown function a(x) in ( 35 ). Once these three functions are determined, all of the pertur- 
bation quantities can be calculated by the appropriate formulae given in this section. 
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5 Numerical results 


In this section, we shall outline a numerical scheme which we have used to integrate the above 
free boundary problem and discuss our numerical results. 

For convenience, we shall drop the hat notation and subscripts c o’ in ( 53 )— ( 57 ). Thus our 
free boundary problem is to integrate 


(l + m) d ( VT dT\ 5T 0T _ 

<7 drf yT + m dr) J + ^ drj dx 


( 58 ) 


subject to the boundary and matching conditions ( 55 )— ( 57 ). The interval (0,oo) is divided 
into three sub-intervals: 


Ti : ( 0 , 771), r 2 : (771, 772), r 3 : (772, 00). 

Our aim now is to integrate ( 58 ) in the intervals and T 3, and iterate on the values of 
771(2), 7 te( x ) ai *d a ( x ) at a gi yen x 80 that the matching conditions ( 56 ) at 771 and 772 are 
satisfied and the integral ( 57 ) over T2 vanishes. 

For the purpose of numerical calculation, it is necessary to work with fixed boundaries so 
in r 1 and T3 we make the transformations 

77 = T7 1 (i)e^, 77 = 772(1)^, ( 59 ) 
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respectively, so that the intervals T i and T 3 now become 


r/:(-oo, 0 ), r 3 ^:(l,oo). 


( 60 ) 


The additional exponentional stretching in (59a) is introduced to accomodate the rapid change 
of T near 77 = 0 (as indicated by (55)). 

In terms of the new variables <f> and ip, (58) becomes 


AT d 2 T dT 

2z— + Ai(x,4>,T)-q^ = Bi(x y <f>,T , —), 


dx 


d<y 


and 


respectively, where 


2x ^ + ^x,^T ) 0 = 5 3 (x ) ^T,^ ) , 


(61) 

(62) 


A 1 = - 

Bi=(l + 


1 + m VT e~ 2 * 


cr T + m rfi 
2xr}[ 1 + m sfT e~ 2 ^ 

nl) 


Vi 


T + m 


dT 

d<t> 


+ 


l + m e 2< ^ m — T ,c?T , 2 
cr r}\ 2 VT(T + m ) 2 d<p 

l + m y/T 1 
a T + m 


(63) 


A 3 = -- 


„ / 217?' \ ,dT 

M 1 + ir ) -* JS9 


dT 1 + 771 1 


m — T 


(— ) 2 

V A ./. ' 


<t 7)1 2 VT(T+m) 2 dip 

Equations (61) and (62) are parabolic partial differential equations, so their solutions can 
be otained by a marching procedure. We shall now use the solution of (61) as an illustrative 
example to explain our numerical scheme. If the values of T, 77 i(x), 772 ( 2 :) and a(x) are known 
at x = i, then the following schem is used to determine these functions at x + e: 

„ Ti-Ti , , T{+i — 2 Ti + Ti_i d ji *n (ka\ 

2x • — f- i4i(x,^j, Ti) • ^2 — Bi\Xj<Pu T l} ^ ), (o4) 


h 2 


where h is the vertical grid spacing, a tilde denotes a quantity evaluated at the position x + e 
and a subscript signifies evaluation at the indicated vertical grid point. In the expression for 
B\y Vi(^) ls replaced by 


ui(s) = 


rh - Vi(x) 


(65) 


where 771 is a guess for 771 ( 2 ; + e). If we replace —00 by <po and use n mesh points in the 77 
direction, we have 

(pi ~ <Pq + ih y <p n = (po + nh = 0 . 
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Application of (64) to i = 1, 2, • - • , n — 1 gives a triagonal matrix equation which can be solved 
after the following boundary conditions are incorporated: 

w 9(1 + m) 1 _ >r_( „ \ (aa\ 


I JL 

fj* . e 4 <h ’ 

1? 

s£ 

11 

(66) 

. The derivative of T at 77 = 771 (x) is 

calculated from 

1 dT 

rfi 8 <f> * =0 — 

^n+ 1 “ T n -i 
rfih 

(67) 

dT 

df 0 

(68) 

- 

■ dr) 


and we define a function f\ by 

/i(d, ffi) — I rj=»7i » (68) 

where the second term is calculated from (34) and a is a guess for a(x + c), which is used in 

the calculation of To(x,?7i) according to (35). Equation (62) can be solved in a similar fashion, 

which leads to a second function /2 : 

, ~ V dcf. 1 dT dT 0 , fiQ x 

/ 2 (a,^) - ” ( 69 ) 

where Tfc is a guess for x + e). For a given guess (a, 771, 772), a third function / 3 is defined by 
means of 

1 J 1 9 # Po«fo, , 5T 0 n 5T 0 1 J ^ 

«-.*.*) -y 4 ^{^(Ka^ ) + ’'W 2 *a7r’’- (70) 

With the aid of the three-dimensional version of the Newton-Raphson method, our program 
iterate on (6,7^,772) until the three error functions become sufficiently small simultaneously. 

The above procedure shows how to march the values of (a(x), Tft(x), *72(2:)) one step forward 
along the stream- wise direction at a given downstream location. The scheme is complete if the 
initial values of (a(x), 771(1), r^x)) are known at a certain initial position x = xo- Such values 
are provided by the weakly nonlinear theory, as we show below. 

The weakly nonlinear theory established in §4 gives the large X(= (x — x n )/e) structure 
for growing Gortler vortices. The present fully nonlinear theory for x — x n = 0(1) should 
then match in the small (x - x n ) limit with that large X structure. Therefore, the desired 
initial conditions for (a(x), 771 (x), ^(x)) and T are imposed near the neutral position x n and 
are obtained by rewritting the large X solutions of the weakly nonlinear theory in terms of the 
original variables x and 77. 

First of all, the initial value of a(x) can be obtained by using the condition that 0 m 0 given 
by (28) vanishes at 77 = 77* so that To in (35), when evaluated at 77 = 77*, can be replaced by 
T j the unperturbed basic state temperature. To determine 77! and 772, we simply have to write 
f = ±C defined between (32) and (33) in terms of the original variables x and 77 with the aid 
of (9), (10), (25) and (30). The result is 

771,2 = v* ± C\j - Zn), (71) 
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where a and b are given by (23); whilst C is determined by solving (44a). Finally, the initial 


value of T can be written as 

T = T + e 3 / 2 0 mO , 


(72) 


where T is again the unperturbed basic state temperature and the other term is the one 
appearing in (14d). Rewriting the latter in terms of the original variables x and 77 with the 
aid of (9), (10), (25), (39) and (42), we obtain 


3 / 2 <U> = ±A 


3Ti b 3/2 


2x n T 2 


(x-x n ) 3 ' 2 -exp(-^f)-U{ 




(73) 


where r - 

_ 2 ya 77 - 77 

\/& \/x — X n 

and where the ‘4-’ and *— ’ signs should be taken for 0 < 77 < 771 and 772 < 77 < 00 respectively. 

In Fig.2, we have shown the evolution of the mean temperature correction e 3 / 2 6 m 0 given 
by (73) downstream of the neutral position x n = 0.5; whilst in Fig.3 we have shown how the 
growth of e 3 ^ 2 6 m 0 depends upon the neutral position x n . We can see that for a fixed value of 
x — i„, e 3 / 2 0 m o decreases drastically with increasing x n . For x — x n = 0.001, our numerical 
calculation shows that e 3 / 2 9 m0 becomes as small as of order 10 -7 when x n = 20. In our 
numerical experimentation, we find that if we choose too large a value for x n , the amplitude 
of the initial Gortler vortex would be too small to have any effect on the evolution of the 
temperature, and as a result, the numerical values of 771 ( 1 ) and 772 ( 1 ) would coalesce into a 
single value as we march downstream. This is why we choose rather small values for x n in 
our following numerical discussion. Such an experimentation also provides a check on our 
numerical scheme as we expect that the two free boundaries would coalesce if no vortices were 
present. 

We now discuss our numerical results. It was found that the above numerical scheme 
converged for sufficiently small values of e and that hi = 0.005, h 2 = 0.004, l - 0.0001 gave 
a stable scheme for the cases investigated and yielded values for 77 ! , 772 and the other flow 
quantities correct to two decimal places, where hi and hi are the vertical grid spacing in rf 
and Tlf, respectively. All cases correspond to a = 0.72, m = 0.509, AT = 1 and to a thermally 
insulated wall for which the basic state solutions have been given in Fu, Hall and Blackaby 
(1990). The first case we considered has the curvature distribution and the neutral position 

given by 

k{x) = (2x) 3 / 2 , x n = 0.4. (74) 

The correponding Gortler number is from ( 1 ) given by G = 21.8044/e 4 . The weakly nonlinear 
theory results (71), (72) and ( 73 ) were used to calculate the initial values of 771 ( 1 ), 772 ( 1 ) 
and a(x) and the initial profile of T at x = 0.401. The numerical scheme described in this 
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section was then used to advance the solution beyond x = 0.401. The numerical values thus 
obtained for 771 (x) and rfe(x) are shown in Fig. 4 and the total temperature distributions at 
x = 0.4, 0.8, 1.2 are shown in Fig.5 with that for x = 0.4 correponding to the unperturbed 
basic state temperature. In Fig. 4 we have also shown the values of ^(x) and 772(2;) calculated 
according to the weakly nonlinear theory result (71), which is strictly valid only for x —x n <C 1. 
In Fig. 6 we have shown the eigenfunction Vq , calculated from (36), at the downstream locations 
x = 0.5, 0.6, 0.7, 0.8. 

The second case we considered corresponds to 

k ( x ) = V2x y x n — 0.3. (75) 

The Gortler number for this case is from (1) given by G = 18.9681/e 4 . This is the case which 
admits a similarity solution in the context of incompressible flows, as has been shown by 
Hall and Lakin (1985). Such a similarity solution is no longer possible here because of the 
contribution to the Gortler number expansion from the basic state curvature, as can be seen 
from (27) and (35). However it can be deduced from (27) that 

H(x) ~ /t(x) = y/2x y as x -4 00, 

and further from (35) that rfr and 772 become independent of x when x becomes large, so that 
a similarity solution is possible for large x. This is verified by our numercal results shown 
in Fig.7 which clearly shows the increasing independence of 77! and 772 on x when the latter 
becomes large. 

6 Secondary instability 

After the large Gortler vortex structure discussed in the previous two sections has been estab- 
lished, we expect that the boundary layer would become susceptible to secondary instability 
of the wavy vortex or vorticity mode type. Thus following Hall and Seddougui (1989), we 
now study the secondary instability of the steady structure described above by superimposing 
spanwise periodic travelling waves on the flow in the two transition layers. We shall confine 
our attention on the upper transition layer; the lower transition can be studied similarly. The 
steady vortex structure in Ha is now our basic state. We signify it by a subscript B and rewrite 
it here for easy reference: 

ub — u + - \-U, vq = v + V y wb ~ W, 

Mi 

PB = P+ -^(P+ -P), T B = f + T, (1) 
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where 


u 


dj_ 

dr,’ 


v = 


/ = 


1 

\PIx 

v- 

T = 


(-Tf + y^I{T )) + \/ 2 ^ 

P . /(g>*7) ■ 

M J /2 M x ’ 

T(x,rj) + o(M°), 


9/(T) 

’ 


( 2 ) 


^ = / o ( z ,*? 2 ) + ^(®, 7 ? 2)-^ 2/3 + c 4 / 3 / 2(*,0 + '"> 

T = fo(x, Tfc) + ^(x, 7ft) • e£ 2 / 3 + 6 4/3 T 2 (x,0 + • . . . 

Here /„ and f 2 are respectively the same as those appearing in (18a) and (38f) and I(T) is 
defined by (5). The harmonic part (U,V,W,P,T) expands as in (38). But without loss of 
generality we rewrite them as 


U = e 4/3 cos ~(U 01 + e 2/3 U n + ••■) + •”. 

€ 

V = e -2 ^ 3 cos — ( Voi + e 2 ^ 3 Vii + ••■) + ’"> 
e 

W = e -1/3 sin ~(W 0 1 + £ 2/3 Wn + •••) + ■’•» ( 3 ) 

e 

T = e 4/3 cos -(Toi + e 2/3 T u + •••) + ‘ * ’ > 

£ 

P = £" 4/3 cos -(Poi + e 2/3 Pu + •••) + •■•» 

€ 

Comparing (3) with (38) shows that U, V, W, T, P here are in turn equal to 2U,2V,2iW,2T,2P 
there. Therefore, the equation satisfied by V 0 i here can be obtained by replacing V 0 i in (48) 
by Voi/2. Thus we have 

+ 5l( I )€^01 = 3^2 ^0 3 1 + ^2( Z )^01> 

We now look for travelling wave solutions superimposed on the above steady state. The 
total flow is then written as 


u = u B + -^-SU*, v = v B + 6V*, w=w B + SW\ 

M\ 

p = p B + SP\ T = T B + ST*, (5) 


where S is a small parameter introduced to facilitate linearization. The linearized perturbation 
equations are obtained by replacing (U,V,W,P,T) in (11)— (15) by ( U + SU*,V + 6V ,W + 
SW*,P + 6P*,T + ST*), putting back the d/dt terms in the momentum and energy balance 
equations, and then linearizing in terms of 6. They are given by 


1_ dU* f dU* 
T dx + Mi dx 


V_9U * 
2x dr) 


- pxr» 


2xTdT) K T dr] 1 split 2 
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+ T^ v '^ + v W )m + WU:] 

t .bv n av J" .... raa , ai/ , J" 

r= l a* 2i a^ vs? 1 r* l ai 21 a, '*' v /5if 1 

— i^.—fB-T^—) - -(£t*—)-u(TU*) - u(T*U ) 

2xTd^T dr, > 2xTdr, K T dr, } M z)z M t)z 

dp vf" 2 tt' i_ a /g \ 1 _ r r 

' /}.» 0~ ’ T3 0~'T a~, I T J'V'T T 2 dt U> 


2x ' T 3 2xT5?7 


1 , 9Vp PjVp_ _ ^dJT _ _ yT _ 4 d n dVp 

T K dx M l dx 2x dr, ’ M “ 2 xT 2 3 ■ {2x)T dr^T dr} } 


1 dP * H{x) T* 2 jW dT * 2£T' 

V^fdr, T 2 e 4 3xT 2 577 + 3v^r 


-^kfk w:) + ? (^ [v ‘ r - + FV -'] + + w ') 

t f dv* r] dv 77 f \ t* fdv t) dv r)f ' N 

T 2 \ dx 2x07? 2xT K j T 2 \dx 2x dr) 2 xT ) 

( 4 dfidv 4A 0,1 WY A fflr, 2 d_l 1 

13 -( 2 x)T 2 dr, dr, T 3 ( 2 x)T 3^25? dr, 3^T dr, z ^ ^ zz ] 

r 4 5A5F* 4A 3 iar A dw ; 

1 3 • (2x)T 2 $77 £77 + 3{2x)TdpT dr] ) + 3 V2^f~d^~ 

1_^VT _l_ ny* \ t — - V ^ -i-^1 

3\/2xf ^ 77 ^ + dr] 3 fXWz ) V2ZT dr, 

\3\/2xT dr, 3 mW *J v^T ^ (^T ^77 + *7 ^ 

+5H - _ J,A (aArr* > ! - + ^ , r 

T 3 e 4 6 x ^77 J \y/2xT dr, + * ) * z 

1 TdV*_ 

+ 'T T 2 ' dt ~ ’ 

1 aw* paw T]_dwp a tp_ A dv; 2^ 

T K dx ^ Ml dx 2x dr, ' spat z 3y/2xT dr, z + 3 • ( 2 x)f 1 


+ V* + _LA ( £™A _ i (_1_ V -9W _J_ y d_W-_ 

*3^ zz + 2xTdr, ( T dr, } T \V2xT dr, + dr, + W ^ z 

+W W) + Z(™L- JL—\ + L (™1 - JL™L\ 

* ' ' j> 2 v dx 2xdr, ) + T 2[ dx 2x dr, > 

, P dV, 1 5A„ , 1 dfidW A d 


+ Ituw + -A_^ + L_^ + -L-^^,_LA ( i^:a T * 

1 3^ zz 3p2xf dr, \/2xT dr, z 2xT 2 dr, dr, 2xT dryT dr, j 
l3 P “ 3^? dr, sfiifdr, z T 2xT 2 dr, dr, ^2xTd77 k T dr, ] ) 
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( 8 ) 


, M_ r _ i_ — 4- v 'I— 4- ** ( i — + f*1 

+ V^T [ V2Zf dr, + z) dr] + s/2if K sfcT dr, z } dr, 

2 dV* 


.dT* 




BW* 


r .s9T 


t w _ _2_£K1 

3 * 3V2xTdr,\ 


vt; + 


4 

-W* 

3 * 3s/2xT dr, 


_ m ,1 T.0W* 

(f T 2 ) Qt ’ 


1 


_/«. + i±-i±vH) + 

K dx Midi 2 xdr,’ K T 2 ’ \/2xT dry T 


v' 


d V* W* 

(ijr) + ~ ^/2xfi 


-<OT v s*" , s«S'-'*r r h" r i*l> 

2v' 


. /_£ ^ ., 2 TT . ^ 1/ _ rr* = o 

+ ^x 2* 3i7 f 3 ^ s/2iT* T 


1 dT* f dT* r, dT* 

T 0x + Mi 0x 


■) + 


2x dr, ’ s/2 xT 

u . uT' #r* 
_ a T “ " 2 xT 2 <7 dr, 


r^* + 


fj£_ 

\2xT 2 


a ( at 


2 xoTdr, K T 


-)jr 


_\7tt / 


+ 


*r * v"+ 


T 

y2 


9T* 
. 3x 


2xTadr, K T dr, 

rwsr + 

r? ar* t' v .\ Wgr ^ar 

7 + //2xf V ) T*\dx 2x 07? s/2xT , 


2x dr, 


1 d^OT* 


' 2xT 2 a dr, dr, 

A 


w-^kl r f-^ (rr * ) - +(rT;w 


a t ar gr ar t,T' 2 tt* 

2xTa^ryT dr, ^ T dr,’ dx 2x T 3 

- 1 _ — ( t£!-TT*] + (i- - Z-)ZE- = 0. 

2^Tdr,\T^ ) +{ T T 2 dt 


(9) 


( 10 ) 


Here /' has been used to denote df/dr, to simplify notation. In the two transition layers, the 
vertical lengthscale is characterized by the variable £ defined by 


( = 


V~V2 

,2/3 1 


( 11 ) 


thus 


JU 0(e - n 


_a 

0f 


= 0(e' 2 ), 


( 12 ) 


where the scales for x and t are deduced from perturbatin equations. Since we are looking for 
travelling wave solutions, it is convenient to define a new variable <f> by 


<j> = t- x 


(13) 
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and use in place of (<p,x,£,z) as the new independent variables. Hence 


d_ d__ d_ d_ d_ 

dx dx dip ’ dt d<p' 


and the operator 


1 + 


£_\d_ 

Mx) dx 


+ 


£ 

dt 


(14) 


which appears in (6)— (10) is transformed to 

dx Mi dcf) 

to leading order. 

Next, we assume that the superimposed disturbance is 7r/2 out of phase with the steady 
state. Therefore, the perturbation quantities take the form 

U* = 6 4 / 3 sin - E- (u 0l + e 2 / 3 u„ + •••) + ••■ + C.C., 

V* = e" 2/3 sin - • E • (v 01 + e 2/3 vn + ■■■) + ••■ + C.C . , 


W* = e~ l ' 3 E ■ (w m0 + e 2 ' 3 «w + •••) + •■• 


4" e cos — • E ■ (twoi + c 2 ^ 3 ttfn + ■■•) + ••• + C.C., 


(15) 


P* = e' 4 ' 3 sin -E- (p* 01 + e 2 ' 3 p*n + ••■) + *•• + C.C., 


T * = e 4 / 3 sin — ■ E ■ (6* 01 + e 2 > 3 6\ x + •••) + ••• + C.C.. 


Here 


E = exp (? f K{x)dx - ^^) ’ (16) 

where fi is the constant frequency and the wavenumber K expands as 

K = K 0 + e 2 ' 3 K 1 +---. (17) 

The scale for <f> in (16) is chosen so that in (14) 

— — = 0 (— ). 

Mi dip K dx J 

We now substitute (15) and (3) into (6)— (10). By equating the coefficients of sin {z/e) 

and €~ 2 Esin(z/e) in (6) we obtain 


7 tT{K o + /oft)^01 + fiQVoi = Oj 

^r(Ko + fott)vn + + foW ~ (*o + 

i 0 ^0 ±0 


( 18 ) 
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( 19 ) 


+ ^0^X1 + TqOIoVqi — 


4/xo d 2 v 0 i 1 dpi 


+ 


01 


3 • (2x)Tq d £ 2 VteTo d( 


H ( x )/a* 2T 0 to* \ I Ao dwo1 1 V- _n 

Equating the coefficients of e~ 7 ^ 3 E cos(z/e) y Ee ~ 7 ^ 3 , and Ee -5 / 3 in (8) gives 

- £(*0 + *<»)«« + ^=f^r ~ pS. - = 0, 


- 7?r{Ko + fotyWmO = 0, 

J-O 


( 20 ) 

( 21 ) 


rj r_^ V — « * 

, Po d 2 w m0 


-M K o + /o«)«»mi -ir \ K 1 + - (JT 0 + /o«)Se| u / m0 


1 , dWoi . TJ . 5u>oi 


2*JJ^ ' r # ) = 0 ' 


( 22 ) 


Finally, from equatig the coefficients of € 4 / 3 E sin( 2 / c) in (9), c 2 / 3 E sin(z/ e) and e°E sin(z/e) 
in (10), we have 

1 m.. 

(23) 

(24) 


1 ^Ol _ n 

71,01 v^fo _ ’ 

J^( K ° + / 0 W 01 + ^/2^2 Vo1 + “£°^01 = °» 

i-(K 0 + /ofiWx + + /o^ - (K 0 + fon&tKi 

1 0 J-0 J o 


+ 


,df 2 2 Tq 

V 2 ZT 2 [ oVl1 + ( d( To ° V °\ 


Ao 

2i<tT^ d£ 2 


H — (Mil + PoT& 6* 01 ) - ^-w m0 9oi — 0. 

< 7 1 0 


(25) 


We now proceed to solve this hierarchy of equations. First, it can be seen that (18) and (24) 
have non- trivial solutions only if 

K 0 = -/ofl. (26) 


It then follows from (16) and (2a) that the travelling wave propagates downstream with 
the same speed to leading order as the that of the basic steady flow. We also note that (21) is 
now automatically satisfied. 

From (24), (23) and (20) we have 


% = -7S&, f 01 ' “’ 01 


1 dv 0i , _ fJo dv 01 
\/2i7o d Z ’ P01 ” y/2xTo d£ ' 


(27) 


With the use of these relations, equation (19) and (25) can be reduced to a 2 x 2 matrix 
equation of the form AC = f for £ = (vn, 0n) T . The inhomogeneous term f involves u 0 x and 
its derivatives. It also involves dti/d^, the expression of which can be obtained from (47) by 
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replacing Voi there by Voi/2. Since the coefficient matrix A has zero determinant, ( has a 
non-trivial solution only if the inner product of f with the left eigenvector of A vanishes. The 
condition then leads, after some algebra, to the equation 

isf - + 3{h + 5 )J ° £ + + / ° fn) 


+ 2jtrs v °' + *“ + 


(1 + q)(2x)7o 


TUmoVoi = 0 . 


The equation involves w m0 as well as Voi- To determine w mQ , we turn to (22) which can be 
shown to reduce to 

d 2 w m o . 2 x% 7>fn\ i 1 ( d2y oi v _ n / on s 

d? l ' fio ^ Kl + fo ^ Wm0 + 2jioTb V d t 2 Vm d ? / °' 9 

With the use of the definitions (49) and (50), equation (28) can be rewriten as 

+ Si(x)£v 0 i = - V^uoi + 52 (x)v 0 i 

+ (1 + ^_ (2 1)70 {{iK^ + iQfU)voi ~ w m0 V 01 } . (30) 

3/x 0 

Finally, after (4) and (30) have been used to eliminate the second order derivatives of Voi and 
voi, equation (29) becomes 

_ ^{iK x + iSlfti)w m0 = {(<JTl + MftiWoiVo! - WnoV 2 ,} . (31) 

Mo °Mo L 

Equations (30) and (31) are to be solved simultaneously to determine the second or- 
der correction K\ to the wavenumber and the frequency fl, subject to the conditions that 
v 01) w m o — ► 0 as 7 ] — ► ±oo so that the travelling waves are confined within the transition 
layers. It is possible to simplify these two equations by scaling the flow properties of the basic 
steady state out of this eigenvalue problem. This can be achieved by introducing the following 
new independent and dependent variables: 


a [2 S l x /3 

v= ~rr^yx'~f^' vou w=Wmo ■ 


In terms of these new variables, equations (4), (30) and (31) become 


-<v = v\ 


^ - (1 + |«n)C« - \i Kv - vV 2 + 2 -sfcVw = 0, 
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(35) 


d? 


V ) 


where 


d( 2 1 + o’ 

ft _ (1 + 


i(Cl( + K)w + 


1 + a 


V 2 w - 


1 4- o_ 

\[§o 2 


i(fiC + K)Vv = 0, 


■Cl, K = 


(1 + a )Tpx _ ( 
PoSl /3 {x) V 


K ' + ^7T' n 

Si( x ) , 


(36) 


Mo'S'iC 1 ) 

Equations (33) -(35) are of the same form as their counterparts for incompressible flows as 
discussed by Hall and Seddougui (1989) and Seddougui and Bassom (1990). In fact, when 
<t = 1, the two eigenvalue problems are identical. The neutrally stable solutions correspond 
to real values of f l and K. Such solutions were first given by Hall and Seddougui (1989) and 
were later improved upon by Bassom and Seddougui (1990). The latter authours numerical 
solution shows that the lowest neutrally stable wavy mode has its eigenvalue pair given by 


(K,Cl) = (0.690,0.372). (37) 

Since Cl = 0 corresponds to a stable wavy mode, when Cl is increased the mode described by (37) 
is more dangerous than any other higher mode because it will occur first. It was conjectured 
by the above authors that in general there will be an infinite number of such neutral solutions. 
This conjecture was further supported by Bassom and Seddougui’ (1990) asymptotic analysis 
which shows that there is indeed a family of neutral modes for X > 1, fi > 1. 

Once the numerical values of K and Cl have been found, the second order correction K x to 
the wavenumber and the frequency Cl can be determined from (36). Note that K\ is a function 
of i. Thus for a given frequency Cl, Ki(x) is the wavenumber for the travelling wave to be 
neutrally stable at x. If Cl is held fixed to be the neutral value at x = x, then ( £ ) will be 
complex when x / x, implying that the travelling wave will experience spatial amplification 
or decay away from the neutral position. 

The numercal values of Ki and Cl are also dependent on the properties of the underlying 
steady state, the solution of which has been shown in sections 4 and 5. We shall not give 
any definite values for K\ and fi for any specific conditions, since the principal aim of the 
present section is to show that neutrally stable travelling wave solutions do exist in the present 
hypersonic context. If necessity comes, the values of K\ and Cl for any specific situation are 
obtainable by using the relevant formulae given in the present paper. 


7 Conclusion 

In this paper, we have given an asymptotic description of the nonlinear development of large 
amplitude Gortler vortices downstream of the neutral positon. We have shown how an asymp- 
totic state can be established under the combined effects of viscosity and nonlinearity. We have 
also investigated the possibility of such a large amplitude vortex structure losing stability to 
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travelling waves of the wavy type. Such an analysis has important applications, for example, 
to the flow in engine inlets and near the control surface of hypersonic vehicles. 

The basis of our present studies is the linear theory given in our previous paper Fu, Hall 
and Blackaby (1990). It has been shown there that taking the large Mach number limit 
has two implications. Firstly, the boundary layer splits into two sublayers: a wall layer and a 
temperature adjustment layer. It is the latterlayer that is most susceptible to Gortler vortices. 
We note that in their studies on the Rayleigh instability, Hall and Cowley (1990), Smith and 
Brown (1990) and Blackaby, Cowley and Hall (1990) found that the temperature adjustment 
layer is also most susceptible to the inviscid instability. Secondly, the boundary layer growth 
has two scales: a short scale related to the similarity variable 77 and the usual scale based 
on the streamwise variable x . The short scale is felt mainly through the 0(M 3 / 2 /(2x) 3 / 2 ) 
curvature of the basic state. Thus in the special case when the wall curvature is proportional 
to ( 2 x)“ 3 / 2 , it exactly counterbalances the basic state curvature and Gortler vortices evolve 
downstream in the same manner as those in incompressible flows. In the more general curvature 
case, boundary layer growth strongly affects the evolution of Gortler vortices and it becomes 
negligible only when the local wavenumber is of 0(M 3 / 8 ) or larger. 

In the present paper, we have confined our attention to the 0(M 3 ^) wavenumber regime 
and thus we have been able to exclude the effects of boundary growth. The neutrally stable 
position is then uniquely defined. The linear theory tells us that when a certain parameter 
is positive, Gortler vortices will grow as they evolve downstream of the neutral position x n . 
In the weakly nonlinear theory presented in section 3, we have determined the evolutionary 
behaviour of growing Gortler vortices in a small neighbourhood of the neutral position where 
Gortler vortices grow at a scale dictated by the variable X = (x — x n )/e. It is shown that the 
mean temperature 0 m o and the first fundamental Vo satisfy two coupled evolution equations; 
whilst the mean streamwise velocity u m0 can be determined from another evolution equation 
once 0 mO and Vo have been found (other first fundamental components are related to Vo). In 
the limit X — * 00 , we have 0 m 0 ~ X 3 / 2 , u m0 ~ X 3 / 2 and Vo ~ X 1 / 2 so that when X = <9(e -1 ), 
that is when x — x n — 0 ( 1 ), the mean temperature and streamwise velocity corrections become 
as large as the basic state. When this happens, the weakly nonlinear theory becomes invalid 
and the further downstream development of Gortler vortices is described by the fully nonlinear 
theory given in section 4. There it is shown that Gortler vortices spread into a region of 0(1) 
depth which is bounded by two transition layers. In the region of vortex activity, the mean 
temperature is determined from a solvability condition for the first fundamentals and thus it 
adjusts itself so as to make any modes neutrally stable everywhere simultaneously. The fact 
that the basic state is now completely altered by the the large amplitude Gortler vortices can 
be seen from (32) and ( 33 ) which show the first fundamental Vq as a forcing to the “modified” 
basic state equations. In the two transition layers viscous effects make the fundamentals decay 
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to zero exponentially, so that above the upper transition layer and below the lower transition 
layer there is only the mean flow. The centres of the two transition layers are determined 
by a free boundary problem, which has been solved numerically in section 5 for a number of 
curvature cases. Thus solutions for the first fundamentals and the mean flow quantities have 
been determined for 0 < 77 < 00 in closed form. 

Once the large amplitude vortex structure described by the fully nonlinear theory has been 
established, transition can be reached by two possible routes in the form of secondary instabil- 
ities, as was shown by Swearingen and Blackwelder (1987). The first secondary instability is 
described here in section 6 which takes the form of time dependent travelling waves confined to 
the two transition layers and which leads to the wavy vortex boundaries observed experimen- 
tally. It is shown that such wavy type secondary instabilities may indeed exist in the present 
hypersonic context. The second possible secondary instability is associated with a Rayleigh 
instability. Relevant results will be given in our next paper. 
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Figure 1. The different flow regions beyond the neutral stability position x n . 
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Mean temperature correction 


ncreasinq 


Figure 2. The growth of the mean temperature correction downstream of the neutral stability 
position x n = 0.4 over a wall with curvature k = (2x) 3 ^ 2 , as predicted by the weakly nonlinear 
theory (5.56). The profiles shown correpond to x=0.6, 0.8, 1.0 and 1.2. 
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Figure 3. The dependence on x n of the amplitude of the mean temperature correction at a 
fixed distance of 0.001 downstream of the neutral stability position x n . The profiles shown 
correspond to x n = 0.4, 0.6, 0.8, 1.0 and to the same curvature distribution as in Figure 2. 
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Figure 4. The development of rfi and with x for the case /c(x) = (2x) 3 / 2 , x n 
— fully nonlinear theory; weakly nonlinear theory. 
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Figure 5 The non-harmonic part of the temperature at different downstream locations for the 
same case as in Figure 4. 
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Figure 7 . The development of rji and 172 with x for the case k(z) = y/2x, 
fully nonlinear theory; weakly nonlinear theory. 
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